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Abstract 

An analytic parametrization of the hyperonic matter equation of state based on microscopic Brueckner- 
Hartree-Fock calculations has been constructed using the realistic Argonne V18 nucleon-nucleon potential 
plus a three-body force of Urbana type, and three models of the hyperon-nucleon interaction: the Nijmegen 
soft-core models NSC89 and NSC97e, and the most recent meson-exchange potential of the Jiilich group. 
The construction of this parametrization is based on a simple phase-space analysis and reproduces with good 
accuracy the results of the microscopic calculations with a small number of parameters. This parametrization 
allows for rapid calculations that accurately mimic the microscopic results, being therefore, very useful from 
a practical point of view. 

PACS numbers: 13.75.Ev, 21.65.Mn, 26.60.-c 
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Neutron stars offer an interesting interplay between nuclear processes and astrophysical observ- 
ables . Properties of neutron stars, such as the mass range, the mass-radius relationship, the 
moment of inertia, the crust thickness or the cooling rate, are closely related to the underlying 
nuclear matter equation of state (EoS) for a wide range of densities and temperatures [jj. Thus, 
its determination is an essential ingredient for understanding such properties. 

At densities near the saturation density of nuclear matter, neutron star matter is thought to 
be mainly composed of neutrons, protons and leptons (electrons and muons) in /3-equilibrium. As 
density increases, new hadronic degrees of freedom may appear in addition to nucleons. Hyperons, 
baryons with a strangeness content, are an example of these degrees of freedom. Contrary to terres- 
tial conditions, where hyperons are unstable and decay into nucleons through the weak interaction, 
the equilibrium conditions in neutron stars can make the inverse process, i.e., the conversion of 
nucleons into hyperons, happen, so the formation of hyperons becomes energetically favorable. 
Although hyperonic matter is an idealized physical system, the theoretical determination of the 
corresponding EoS is an essential step towards the understanding of properties of neutron stars. 
Moreover, the comparison of theoretical predictions for the properties of these objects with the 
observations can provide strong constraints on the interactions among their constituents. Since the 



pionneering work of Ambartsumyan and Saakyan [4] the EoS of hyperonic matter has been con- 
sidered by several authors either from phenomenological js-12] or microscopic 13-18] approaches. 

In phenomenological approaches the input is a density-dependent interaction which contains 
a certain number of p arameters adjusted to reproduce experimental data. Within this approach 
Balberg and Gal Hj derived an analytic effective EoS using density-dependent baryon-baryon 



potentials based on Skyrme-type forces including hyperonic degrees of freedom. The features of 
this EoS rely on the properties of nuclei for the nucleon-nucleon (NN) interaction, and mainly 
on the experimental data from hypernuclei for the hyperon-nucleon (YN) and hyperon-hyperon 
(YY) interactions. This EoS reproduces characteristic properties of high-density matter found in 
theoretical microscopic models. Within the same scheme, several authors 19|, [20j have developed 
Skyrme-like YN potentials to study properties of single- and multi-A hypernuclei with the Skyrme- 
Hartree-Fock formalism. 



An alternative phenomenological approach invo. 
mean field theory (RMFT) of interacting hadrons 



ves the formulation of an effective relativistic 



2l| . This fully relativistic approach treats the 



baryonic and mesonic degrees of freedom explicitely, and is, in general, easier to handle because it 
only involves local densities and fields. The EoS of dense matter with hyperons was first described 

n nn 

within the RMFT by Glendenning [5] and latter by other authors [6H10I]. The parameters in this 
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approach are fixed by the properties of nuclei and nuclear bulk matter for the nucleonic sector, 
whereas the coupling constants of the hyperons are fixed by symmetry relations, hypernuclear 
observables and compact star properties. 

In microscopic approaches, on the other hand, the input are two-body baryon-baryon interac- 
tions that describe the scattering observables in free space. These realistic interactions have been 
mainly constructed within the framework of a meson-exchange theory, although recently a new 
approach based on chiral perturbation theory has emerged as a powerful tool. In order to obtain 
the EoS one has to solve the complicated nuclear many-body problem 
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23J. A great difficulty 



of this problem lies in the treatment of the repulsive core, which dominates the short-range be- 
havior of the interaction. Various methods have been considered to solve the nuclear many-body 



problem: the variational approach 



241 ]. the correlated basis function (CBF) formalism [25f] , the 



self-consistent Green's function (SCGF) technique 26], or the Brueckner-Bethe-Goldstone (BBG) 
[27 1 and the Dirac-Brueckner-Hartree-Fock (DBHF) theories 28]. Nevertheless, although all of 



them have been extensively applied to the study of nuclear matter, up to our know 
BBG theory 



eel 



ge, only the 



la], have been 



13H16I] , and very recently the DBHF one [17| , and the V\owk approach 
extended to the hyperonic sector. 

The microscopic approach is in general technically complex and very time consuming. Therefore, 
from a practical point of view, it would be interesting to have an analytic parametrization of the 
hyperonic matter EoS based on such approach that allow to mimic the microscopic results in a fast 
way with a small number of parameters. In the present work we will build a density functional for 
the EoS based on microscopic Brueckner-Hartree-Fock (BHF) calculations of hyperonic matter. 
In addition to the nucleonic degress of freedom we will consider only A and Yr hyperons in the 
construction of our functional, the reason being that these two types of hyperons are the ones 
app earing first in calculations of /3-stable neutron star matter based on microscopic approaches 
16]. The other hyperons, E°, E + , H° and S~, being heavier, either do not appear or only show 



BE 



up at very large densities in microscopic calculations. 

After a brief review of the BHF approach to the EoS, we will detail the construction of the 
functional. We will finish by testing and discussing the quality and validity of our parametrization. 

Our calculations are based on the BHF approximation of the BBG theory extended to the 
hyperonic matter case 



13Ml6|] . Therefore, our-many body scheme starts by constructing all the 
baryon-baryon G matrices, which describe in an effective way the interaction between two baryons 
in the presence of a surrounding medium. The G matrices can be obtained by solving the Bethe- 
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Goldstone equation, written schematically as 

EQBiB 
B . B , uj-E Bi - E B . + ir) 

where the first (last) two subindices indicate the initial (final) two-baryon states compatible with 
a given value S of the strangeness, namely NN for 5 = and YN for S = — 1; V is the bare 
baryon-baryon interaction (NN or YN); Q BiBj is the Pauli operator, that prevents the intermediate 
baryons Bi and Bj from being scattered to states below their respective Fermi momenta; and the 
starting energy u: corresponds to the sum of the nonrelativistic single-particle energies of the 
interacting baryons. We note here that, although we have considered only A and X - hyperons in 
the construction of our parametrization of the EoS, S° and S + hyperons have been also taken into 
account in the intermediate YN states when solving the Bethe-Goldstone equation. The interested 
reader is referred to Refs. 13141a] for computational details. 



The single-particle energy of a baryon Bi is given by 

E Bi {k) = M Bi + ^- + U Bi {k). (2) 

Here M Bi denotes the rest mass of the baryon, and the single-particle potential U Bi represents the 
averaged field "felt" by the baryon owing to its interaction with the other baryons of the medium. 
In the BHF approximation, U Bi is given by 



U Bi $) = Y,U { ^\k) = Re ^^n B ,(|P|)(fcP|GH BiS , jBiB ,(o; = E Bi {k) + E Bj {k'))\kk') , (3) 

b j k 



where a sum over all the different partial contributions, U B 3 (k), is performed, n Bj (\k\) is the 
occupation number of the species Bj, and the matrix elements are properly antisymmetrized when 
baryons Bi and Bj belong to the same isomultiplet. We note here that the so-called continuous 
prescription has been adopted for the single-particle potentials when solving the Bethe-Goldstone 



equation, since, as shown by the authors of Refs. 



29 



30( | . the contribution to the energy per particle 



from three-body clusters is minimized in this prescription. 

All the calculations carried out in this work have been performed with the realistic Argonne 
V18 [31] NN interaction supplemented with a NNN three-body force (TBF) of Urbana type which, 
for use in BHF calculations, was reduced to a two-body density-dependent force by averaging 



over the third nucleon in the medium 



32[ | . This TBF contains two parameters that are fixed by 



requiring that the BHF calculation reproduces the energy and saturation density of symmetric 
nuclear matter. The interested reader is referred to the works of Refs. 33l43a| for a recent analysis 
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of the use of TBF's in nuclear and neutron matter. The YN G-matrices have been constructed 
using three YN models: the Nijmegen soft-core models NSC89 3(| and NSC97e [3?]], and the 
most recent YN meson-exchange potential of the Jiilich group [3a ]. In the following we will use 
the names Jiilich, NSC89 and NSC97 models to denote the three different NN+NNN+YN models 
under consideration, since the pure nucleonic part is the same in all of them. We note that the YY 
interaction has not been considered in the present work due to the large uncertainties still existing 
in this sector. 

Once a self-consistent solution of Eqs. (H|) and ([3]) is obtained, the total energy per particle is 
easily calculated: 



2M B , 2 



T V 



This quantity is a function of the particle densities p n ,p p ,p\, and or, equivalently, of the 
total baryonic density p = p n + p p + p\ + /?£-, the hyperon fraction Y = (p E - + pa)/p, the 
isospin asymmetry (3 = (p n — p p )/(p n + Pp), and the asymmetry between the A and S~ hyperons 

" = (Ps- ~ Pa)/(pe- + PA)- 

As referred before, Brueckner-type calculations are very time consuming since one has to solve 
a self-consistent set of coupled-channel equations for different strangeness sectors. Therefore, from 
a practical point of view, it would be interesting and useful to characterize the dependence of 
the total energy per particle E/A on the particle densities p n ,Pp,PA, and p^-, or, alternatively 
on p, Y, j3 and a, in a simple analytical form. The free Fermi gas contribution, T/A, is already 
analytic, reading 

T ^ 3 fi 2 k 2 Fi Pi _ 3 h 2 k 2 F 1 



A ^ 5 2Mi p 5 2 2 

i=n,p,A,E 



_L (1 _ y)5/3 (1 + 05/3 + _L (1 _ y)5/3 (1 _ ^5/3 
+ ±. yV3 (1 _ a) 5/3 + J_ F 5/3 (1 + a) 5/3 

Ma M s - 



(5) 



where kp i = (S^pi) 1 / 3 and we have defined hp = (37r 2 p/2) 1 / 3 . An idea of the possible terms 
appearing in the correlation energy contribution, V/A, can be obtained from the following phase 



39J for isospin 



space analysis of the single-particle potentials, similar to the ones perfomed in Ref. 
asymmetric matter and in Ref. for spin-polarized isospin asymmetric matter. Replacing the 
matrix elements (kk'lG^BiB^BiB^i^ = E Bi + E Bj )\kk') by an average value, gSiBjik, p,Y, f3,a), 
in the Fermi sphere with radius k! < k F j , and integrating over the corresponding Fermi sea, the 
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(6) 



single-particle potentials of the four baryon species under consideration can be written as 

Un(k) ~ g-nnPn + gnpPp + QnAPK + ffnS" P£" , 

U p (k) ~ 5f pn /) n + gp P /)p + 5 P APA + 5 p s-Ps- i 

^a(^) ~ SAraPra + SApPp , 
^S" 0) ~ 5s-nPn + SE"pPp ■ 

For small values of the hyperon fraction, the isospin asymmetry, and the asymmetry between A's 
and S~'s, one can neglect the dependence on Y, (3 and a of the average G matrices assuming 
gBiBj {k, p,Y,p,a) ~ ffBiBj {k, p) , and 

5nn « 5pp = gi(k,p) , 
Pnp ~ 5pn = 92(k,p) , 

5nA ~ 5An ~ 5pA « SAp = 33 P) , ( 7 ) 
fnS" « 3s-n = 94,{k,P) , 

5pE- ~ fe-p = 05 • 

Note that the quantities <7i(A;, /?) receive contributions from different isospin (T) and strangeness (S 1 ) 
channels. Whereas receives contributions only from the isospin triplet and zero strangeness 

channel, g2(k,p) has in addition a contribution from the isospin singlet, g^(k,p) and g^(k,p) are, 
respectively, purely isospin 1/2 and 3/2 with strangeness —1, and g*,{k,p) has contributions from 
T = 1/2 and T = 3/2 with S = -1. 

Using the set of Eqs. (0) and 0, the single-particle potentials can then be written as 



P 



U n (k) ~ £ (1 - Y) gi(k, p) (1 + (3) + ^(fc, p) (1 - /3) 



Op(*0 ~ ? ( X " y ) S2( fc > p) (1 + P) + gi(k, p) (1 - P) 



+ P -Y 
2 



+ P -Y 
2 



93(k, p) (1 - a) + g±(k, p) (1 + a) 



9z{k, p) (1 - a) + g 5 (k, p) (1 + a) 



(8) 



(9) 



U A (k)~pg 3 (k,p)(l-Y) 



(10) 



^s- (*0 ~ 7; (1 - 3^) p) (1 + /3) + 55 (fc, P) (1 - P) 



(11) 



where the particle densities p n , p p , p\ and p%- have been written in terms of p, (3, Y and a. These 
equations show explicitely the dependence of the single-particle potentials on the hyperon fraction, 
the isospin asymmetry and the asymmetry between A's and E~'s. This dependence is tested in 
Fig. [1] for the Julich model where the value at k = of the single-particle potentials TJ n ,U p ,U\ 
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and E/e- at the saturation density (po = 0.175 fm~ 3 in our model) is plotted as a function of Y 
(with f3 = a = 0), f3 (with Y = a = 0) and a (with Y = 0.1,(3 = 0), in the left, middle and 
right panels, respectively. Similar dependences has been obtained also for the NSC89 and NSC97 
models. The above equations predict a linear variation of the single-particle potentials with respect 
to Y, (3 and a. This prediction is quite well confirmed from the microscopic results reported in Fig. 
[H although deviations from the linear behavior are found at higher values of Y, f3 and a. These 
deviations have to be associated to the dependence of the average G matrices on Y, (3 and a, which 
has been neglected in the present analysis (see the set of Eqs. (J7J)). 

Now, using Eqs. (f8j)- (fTTj) . and replacing the quantities gi(k,p) by their averages cji(p) in the 
corresponding Fermi spheres, one can see, after integration, that the correlation energy behaves 
like 

V ffi(p) ( 2 , 2\ , h(p) , MP) , , \ 9i{p) | MP) /i,! 

T ~ -= (Pn + Pp) + PnPp + [fin + Pp) PA H PnPz~ + PpPY,~ , (12) 

A 2p F/ p p p p 

or, replacing the particle densities in terms of p, f3, Y and a 

\ ~ P -UP) (1 " y? (1 + /? 2 ) + \UP) (1 - Yf (1 - /3 2 ) + P -MP)Y (1 -Y)(l- a) 

+- A HP)Y (1 ~Y)(1 + a) (l + /3) + £g s (p)Y (1 - Y) (1 + a) (1 - /3) . (13) 

From this simple analysis we can finally infer the form of the correlation energy 

j = V l {p) (1 - Y) 2 (1 + /3 2 ) + V 2 (p) (1 - Yf (1 - f) + V 3 (p)Y (1 - Y) (1 - a) 

+y 4 (/>)y (i - y) (i + «) (i + 0) + vup)^ (i - y) (i + a) (i - /?) . (u) 

The coefficients V\ (p) , V2 (p) , V3 (p) , V4 (p) and V5 (p) have been fitted to reproduce the microscopic 
BHF results corresponding to the following five set of values of Y, f3 and a: (Y = 0,j3 = 0, a = 0), 
(Y = 0,0 = 1, a = 0), (y = 0.1,/3 = 0.875, a = 1), (Y = 0.15, /3 = 0.7, a = 0.5), and (Y = 0.2, /3 = 
0.5, a = 0). The first two sets guarantee that the parametrization reproduces the microscopic 
results for symmetric nuclear matter and pure neutron matter. The other three have been chosen 
in order to mimic three representative /3-stable matter compositions for densities above the hyperon 
threshold obtained with microscopic approaches (see e.g., Refs. Q and Q). It is clear that the 
determination of these coefficients is not unique. However, we have checked that with the choice 
of this set of values of Y, (3 and a, we get a parametrization that reproduces with good quality (see 
Fig. 0] and the discussion below) the results of the BHF calculation for a wide range of arbitrary 
values of Y, (3 and a. In addition, we have adjusted the density dependence of the coefficients Vi(p) 
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in the following functional form 

Vi(p) = api + bp s i = l-5 , (15) 

where the set of parameters a, 7, b and 5 for the five coefficients are given in Tables HT lIIIl The 
density dependence of the coefficients Vi(p) together with the functional defined in the above 
equation is shown in Fig. [2] for the three models considered. It is seen that the functional of Eq. 
(|15|) reproduces reasonably well the microscopic results for the three models in the whole range of 
densities explored (0.01 fm~ 3 < p < 0.5 fm~ 3 ). 

In order to test the quality of our parametrization, we show for the three models, in the left, 
middle and right panels of Fig. O the correlation energy at p = po as a function of Y (with 
(3 = a = 0), /3 (with Y = a = 0) and a (with Y = 0.1,(3 = 0). Circles, squares and triangles 
show the result of the microscopic BHF calculation obtained with the Jiilich, NSC89 and NSC97e 
models, respectively. From the figure it can be seen that the dependence on Y, f3 and a predicted 
by Eq. ()14p (quadratic in Y and /3, and linear in a) is well confirmed from the microscopic results. 
For completeness, we finally compare in Fig. [H the results for the correlation energy as a function 
of the density for three arbitrary sets of values of Y, /3 and a: (i) Y = 0.08, f3 = 0.6, a = 0.6, (ii) 
Y = 0.15, /S = 0.2, q = 0.2, and (iii) Y = 0.17, f3 = 0.4, a = 0.75, obtained from the microscopic 
BHF calculation (symbols) and from the parametrization (lines) for the three models considered. 
The quality of the parametrization is quite good for the three models, as it can be seen in the 
figure, with deviations from the microscopic calculation, at the higher density explored, of at most 
2% and 5 — 6% for the NSC97 and Jiilich models, respectively, and of about 10 — 11% for the 
NSC89 one. These deviations, as it has been said before, have to be associated to the dependences 
on Y, f3 and a neglected in the construction of the parametrization. 

To summarize, we have constructed an analytic parametrization of the hyperonic matter equa- 
tion of state based on microscopic Brueckner-Hartree-Fock calculations using the realistic Argonne 
V18 NN potential plus a NNN TBF of Urbana type, and three models of the YN interaction: the 
Nijmegen soft-core models NSC89 and NSC97e, and the most recent meson-exchange potential 
of the Jiilich group. The construction of this parametrization is based on a simple phase-space 
analysis, and reproduces with good accuracy the results of the microscopic calculations, allowing 
for rapid calculations that accurately mimic the microscopic BHF results, being, thus, very useful 
from a practical point of view. Our parametrization will be extended in a future work to include 
finite temperature effects, necessary to describe the properties of newborn neutron stars, and the 
conditions of matter in relativistic heavy-ion collisions. 
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FIG. 1: (Color online) Neutron (circles), proton (squares), A (triangles) and S~ (diamonds) single-particle 
potentials at k = and p = po as a function of Y (left panels), j3 (middle panels) and a (right panels) 
obtained with the Julich model. Upper (lower) panels show the results for neutrons and protons (A and 
£-). 
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FIG. 2: (Color online) Density dependence of the coefficients Vi{p) of Eq. (TH| . Circles, squares and triangles 
show the result of the microscopic BHF calculation obtained with the Julich, NSC89 and NSC97e models, 
respectively, whereas solid, dotted and dashed lines refer to the parametrization defined in Eq. (|15p and 
Tables EHUD 
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FIG. 3: (Color online) Correlation energy V/A at p — po as a function of Y (left panel), j3 (middle panel) and 
a (right panel). Circles, squares and triangles show the result of the microscopic BHF calculation obtained 
with the Julich, NSC89 and NSC97e models, respectively. 
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FIG. 4: (Color online) Correlation energy V/A as a function of the densty for three arbitrary sets of 
values of Y,/3 and a: Y = 0.08, /J = 0.6, a = 0.6 (circles), Y = 0.15, /J = 0.2, a = 0.2 (squares) and 
Y = 0.14, P — 0.4, a = 0.75 (triangles). Symbols show the result of the microscopic BHF calculation, 
whereas dashed lines refer to the parametrization defined according to Eqs. (| 14[) and (fl~5|) . and Tables IT MlIIt 
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Coefficient 


a 7 b S 


Vi 


-65.0189 0.843983 166.944 1.89579 


v 2 


-144.122 0.628802 82.4707 0.829031 


v 3 


-241.211 0.984562 195.95 1.99311 


v 4 


-123.882 0.999992 76.3707 0.90001 


v 5 


-98.8994 0.826319 14.7032 0.830927 



TABLE I: Set of parameters a, 7, b and S characterizing the density dependence of the coefficients Vi(p) for 
the Jiilich model. The parameters 7 and 8 are dimensionless, whereas the units of a and b are MeV x fm 37 
and MeV x fm 3<5 , respectively. 
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Coefficient 


a 7 b 6 


Vi 


-65.0189 0.843983 166.944 1.89579 


v 2 


-144.122 0.628802 82.4707 0.829031 


v 3 


-142.908 0.897287 183.341 1.99791 


v A 


-123.583 0.987746 72.1931 0.983861 


v 5 


-114.735 0.629932 57.809 0.632833 



TABLE II: Same as Table U for the NSC89 model. 
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Coefficient 


a 7 b 6 


Vi 


-65.0189 0.843983 166.944 1.89579 


v 2 


-144.122 0.628802 82.4707 0.829031 


v 3 


-108.951 0.71094 149.79 1.8744 


v 4 


-95.221 0.999994 11.6305 0.901149 


v 5 


-21.5452 0.34935 76.0311 1.88421 



TABLE III: Same as Table U for the NSC97e model. 
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